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ABSTRACT 

A number of metal-rich white dwarfs (WDs) are known to host compact, dense particle disks, which 
arc thought to be responsible for metal pollution of these stars. In many such systems the inner radii 
of disks inferred from their spectra are so close to the WD that particles directly exposed to starlight 
must be heated above 1500 K and are expected to be unstable against sublimation. To reconcile this 
expectation with observations we explore particle sublimation in H-poor debris disks around WDs. 
We show that because of the high metal vapor pressure the characteristic sublimation temperature 
in these disks is 300 — 400 K higher than in their protoplanetary analogues, allowing particles to 
survive at higher temperatures. We then look at the structure of the inner edges of debris disks and 
show that they should generically feature superheated inner rims directly exposed to starlight with 
temperatures reaching 2500 — 3500 K. Particles migrating through the rim towards the WD (and 
rapidly sublimating) shield the disk behind them from strong stellar heating, making the survival 
of solids possible close to the WD. Our model agrees well with observations of WD+disk systems 
provided that disk particles are composed of Si-rich material such as olivine, and have sizes in the 
range ~ 0.03 — 30 cm. 

Subject headings: White dwarfs — Accretion, accretion disks — Protoplanetary disks 



1. INTRODUCTION. 

About two dozen (at the moment of writing) white 
dwarfs (WDs) are known to exhibit near-IR excesses 
in their spectra (e.g. [Zuckerman fc Bccklin 1987*; 
Kihc et all I200I 2006; lJura et al 2 0,09 , Farihj ct aU 
20091). This is us ually interpreted ([Graham et al.l ll990l: 
Jurali2003al 2006; iFarihi et all 120091 ) as evidence for the 
existence of nearby solid debris reprocessing stellar radi- 
ation in the IR. Detailed spectral modeling of excesses 
generally supports the idea that debris particles are ar- 
ranged in a disk-like configuration, which is optically 
thick but geometrically very thin'^, thus h aving prop- 
erties very similar to the rings of Saturn (jCuzzi et al.l 
|2010[) . These disks are relatively compact, < 1 R©, ex- 
cluding the possibility of a primordial origin because of 
the long cooling ages of the WDs around which they 
are fou nd, ~ 0.1 — 1 Gyr. It was prop osed bv | Jura 
( 2003bfl ■ following an earlier suggestion bv lAlcock et al. 
(|1986( ). that these disks are produced by tidal disruption 
of asteroid-like bodies launc hed on low-periastron orbit s 
by distant massive planets (jDebes fc Sigurdssonll2002[ ). 
which have survived the asymptotic giant branch (AGB) 
phase of the evolution of the central star. 

All WDs possessing compact debris disks exhibit at- 
mospheres which are polluted (sometimes heavily) with 
metals, putting these WDs into the DAZ and DBZ classes 
(jFarihil |2011[ ). This observation strongly suggests that 
many (if not all) metal-rich WDs are polluted by accre- 
tion of high-Z elements from the compact debris disks 
around them. In cases where near-IR observations do 
not reveal the presence of a conspicuous disk of solids, 
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the disk may simply be too tenuous to reprocess enough 
stellar radiation to make itself visible. Alternatively, the 
disk of solids may have dispersed some time ago but the 
metals can still be present in t he WD atmosphere be- 
cause of their long settling time ()Metzger et al.|[20T^ . 

In this scenario of WD metal pollution (which is cer- 
tainly valid for systems with known disks) the issue of 
metal transfer onto the WD surface must be addressed, 
as the disk of solids cannot extend inward all the way 
to the stellar surface at Because of the high effec- 
tive temperature of these WDs, ~ (7 - 20) x 10^ K, 
solids must sublimate at some inner radius i?m produc- 
ing metal gas, which is subsequently accreted onto the 
WD via a conventional accretion disk. The existence of 
inner cavities in disks of solid debris follows directly from 
the shape of their spectral energy distributions (SEDs), 
which generally show a lack of emission corresponding to 
temperatures in excess of ~ 200 K. 

It has been shown b y I Rafikovl ()2011af ) and 
IBochkarev fc Rafikovl (|2011l ) that Poynting-Robertson 
(PR)drag on the disk of particles naturally drives 

accretion of solids at rates Mz ^ 10^ — 10® g s~^. 
Their sublimation feeds metal gas accretion onto the 
WD surface at the same rate. Even higher Mz can be 
achieved if the disk of solids can couple v ia aerodynamic 
drag to the surrounding gaseous d isk fCansicke et alJ 
120061 2007; iBrinkworth et all 120091: uMclis ct al. 2010), 
which naturally for ms via sublimation of particles at 
R^n tRafikov..201 1b; IMetzgereFaLllMl . 

1.1. Inner rim puzzle 

Existing debris disk models used to fit SEDs try to link 
the radius of the inner edge of the disk i?i„ to a certain 
value of the "sublimation temperature" Tg, which de- 
pends on physical properties of the constituent particles. 
In these models particles sublimate at radii where their 
temperature T exceeds the sublimation temperature T^, 
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and the inner radius Rin corresponds to T(_Ri„) = T^. 
Then the determination of Rin hinges upon the proper 
choice of Tg and the knowledge of a relation between T 
and r. 

The characteristic value of Tg usually assumed for 
disks around WDs is 1300 - 1500 K. This is a typical 
sublimation temperature for the Si-rich solids, such as 
olivine o r calciu m-aluminum inclusions (CAIs) based on 
iLoddersI (|2003f ) who calculated condensation tempera- 
tures of different species in the proto-Solar disk assum- 
ing Solar abundances of elements. There is good evi- 
dence that Si-rich material indeed represents a significant 
fraction of mass in the debris disks around WDs, both 
from the detectio ns of the Si featu re at 10 /zm in spec- 
tra of such disks ()Jura et alj 12009 ') and the atmospheric 
compositions of t heir ho st WDs (_Zu cker man et alJl2007l : 
iKlein et aUlMol 2011: lJura et al.li2012f ). However, as 
we show in ^ this estimate of Tg needs to be seriously 
revised for the typical conditions in circum-WD debris 
disks. 

There is certain ambiguity regarding the equilibrium 
temperature of the disk particles. In Appendix |^ we 
consider their thermal balance by looking at different 
heating and cooling processes, and find stellar heating 
and radiative cooling of particles to dominate the bal- 
ance. In this case particles directly exposed to starlight 
or located in the optically thin parts of the disk, are 
heated to a temperature 
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In particular, this estimate is appropriate for particles at 
the inner edge of the optically thick disk because these 
are directly illuminated by the star. However, behind 
the narrow rim of directly exposed particles the disk is 
illuminated by starlight only at i ts surface at a g razing 
incidence angle ( « (4/37r)i?i,/r (iFriedjungl 119851 ) for a 
geometrically thin (i.e. flat) disk. This is because of the 
shielding that is provided by the rim particles against di- 
rect starlight for the disk just outside the rim, see Figure 
[1] for illustration. The equilibrium temperature of par- 
ticles in the optically thick parts of the disk is given by 
(IChiang fc Goldreichl[l997l ) 
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This expression is valid in the shielded parts of the disk 
where the near-IR emission is produced. 

The assumption that particles sublimate at a single 
temperature Tg implies that di rectly exposed rim parti- 
cles cannot be hotter than Tg (Raflkovl l2011al , b). How- 
ever, using equation ([1]) to estimate T(r) at the inner 
edge of the disk and taking Tg « 1500 K one finds that 
the near-IR contribution to the SED produced by the 
disk is too weak to account for observations. Indeed, 
combining equations ([1]) and ([2]) one finds 
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Thus, when Tthin is close to the sublimation tempera- 
ture the temperature Tthick in the shielded part of the 



disk must be substantially lower than Tg (since necessar- 
ily Tthin < T^)- For example, taking T^ = 10^ K and 

Tthin (i?m) = T« = 1500 K one finds Tthick (i^^r^) « 660 K. 

As a result, the SED of such a disk is going to be very 
deficient of the near-IR flux corresponding to emission at 
temperatures > 1300 — 1500 K. However, disk SEDs typi- 
cally exhibit considerable emission by material heated in 
excess of 1000 K, see Table [2l This is hard to reconcile 
w ith only the inner rim of the disk being heated to Tg. 

lJura fc Xul ([201^ suggested that this problem can be 
resolved if the inner edge of the disk is set by sublima- 
tion occurring in the optically thick part of the disk illu- 
minated by the star at grazing incidence, rather than in 
the thin inner rim of directly exposed particles. This is 
equivalent to determining the value of Rin by using the 
expression ^ instead of ([T]) in equation T(i?i„) = Tg. 
If that were true, however, then the temperature at the 
inner rim must be higher than T,, see equation ([3]), and 
rim particles would be sublimating, exposing the parti- 
cles behind them to direct starlight. As a result, the 
inner rim would recede to larger distance from the WD 
until it reaches the radius where Tthin = Tg, so we go 
back to the previously considered situation with its in- 
trinsic problems. Only this configuration is going to be 
in stable phase equilibrium as long as sublimation is ide- 
alized as a step-like process, i.e. that particles turn into 
gas as soon as they reach Tg. Such an equilibrium was 
assumed in Rafikov (2011a, b) to determine Rin. 

This set of conflicting arguments suggests that our un- 
derstanding of the location and structure of the inner rim 
is in some ways incomplete. The goal of the present work 
is to fill these gaps and to provide a more detailed picture 
of the sublimation of solids at the inner edge of the disk 
by focussing on two effects. First, in fJD we show that 
sublimation in hydrogen-poor debris disks around WDs 
is different from sublimation in the proto-Solar disk re- 
sulting in Tg being higher than 1500 K. Second, we show 
in f}3] that particle sublimation at the inner rim of an 
optically thick disk is a dynamic process, which makes 
it possible for the rim particles to reach temperatures in 
excess of Tg before sublimating. In 21 we look at subli- 
mation in optically thin disks. We apply our theory to 
observed disk-hosting WDs in ij5l and discuss our find- 
ings in ij6l A summary of our main results can be found 
in gZl 

2. SPECIFICS OF SUBLIMATION IN THE 
CIRCUM-WD DISKS. 

Solid particles of certain composition surrounded by 
vapor with the same elemental abundance grow by con- 
densation of molecules or atoms arriving at their surfaces 
from the gas phase and lose mass due to sublimation. For 
a particle of mass rup and surface area Sp surrounded by 
vapor at pressure Pyap one can write the following mass 
evolution equation (jCuhathakurta fc Drainc 1989i): 
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where /i is the mean molecular weight of the particle ma- 
terial. Here the first term in brackets describes conden- 
sation {{a) is the accommodation coefficient — sticking 
probability of gas particles impacting the solid surface). 



3 



while the second term characterizes subhmation from the 
particle surface (m is the mass loss rate per unit surface 
area due to sublimation). 

When the vapor pressure becomes equal to the satu- 
rated vapor pressure P^^p(T) at a given temperature T, 
the equilibrium between the loss and gain processes is 
established and dnip/dt — 0. This allows us to express 



m(T) = (a)P-*(r) 



2TrkBT 
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The concept of sublimation temperature Ts implies the 
process of phase transition from solid to gas to occur in a 
step-like fashion. At T = Ts the vapor saturates and an 
infinitesimal increase of temperature leads to slow (quasi- 
static) conversion of solid into gas. Then one can again 
assume dmp/dt — and the right-hand side of equation 
(U) then provides us with an implicit relation for Tg as 
a function of the vapor pressure Pvap, as long as the 
dependence rri(T) is known. 

Considerations based on the Clausius-Clapeyron rela- 
tion suggest that 
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where the constants /3 and To are specific to a particu- 
lar particle composition. Since the strongest dependence 
of P^^p on T occurs through the exponential factor, the 
power-law dependence on T in equation ([5]) can be ab- 
sorbed into the (approximately) constant pre-factor, so 
that m{T) is approximated as 

m(T) = (a)i^oe-^°/^, (7) 

where Kq is a constant. This is the form of m{T) that 
we adopt in this work. 

In the following we will consider a number of different 
materials which can represent the composition of disk 
particles. We summarize the parameters Kq and Tq for 
different species considered in this work in Table [TJ and 
provide details of their calculation in Appendix |B] 

Setting the left hand side of equation (|4]) to zero and 
using m in the form ^ we obtain the (implicit) depen- 
dence of Tg on the vapor pressure: 
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with As ^ 1. According to this expression Ts is higher 
for larger Pvap, even though the dependence is rather 
weak (logarithmic). 

The vapor pressure in the gas around the rim can be 
easily estimated if the disks consist of particles with iden- 
tical composition. In general this does not have to be 
true but we still adopt this assumption for simplicity. 
Then metals detected in the WD atmosphere come from 
accretion of this material in the gas phase, and the mea- 
surement of the corresponding mass accretion rate 
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(where v = a^c^/il is the kinematic viscosity, a^, is the 
effective viscosity parameter, D, is the Keplerian angular 
frequency, and Cg and Eg are the sound speed and the 



Table 1 

Sublimation properties of different materials 



Material 


Ko 


To 




Tsub(l dyne cm 2) 




-1 -9 -1 

g cm s 


K 


mp 


K 


Olivine 


1.6 X 10^ 


68100 


141 


2100 


Graphite 


9.2 X 10^ 


81200 


12 


2600 


CAI 


1.1 X IQio 


69400 


274 


2000 


Iron 


2.3 X 10^ 


45400 


56 


1600 


AI2O3 


8 X 10^ 


80500 


102 


2300 


SiC 


6 X 10* 


73700 


40 


2300 



surface density of the gas) provides an estimate of the 
vapor pressure: 
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Here Mz,s = Mz/(10« g s-i), Af.^i = M^/Mq, a,,_2 = 
and Cs,i = Cs/(1 km s^^) is the characteristic 
value of the sound speed for Si-rich material heated to 
temperature of several 10'^ K. According to equation 
Pvap ^ 1 dyne cm~^ for Mz ^ 10* g s~^, which is a 
characteri stic mass a ccretion rate of metals due to the PR 
drag (Rafi kovll2011a[ ). It should also be mentioned that 
a number of metal-rich WDs with debris disks exhibit 
much higher values of Mz, easily reaching 10^ — 10^° g 
s~^. In these systems Pvap ^ 10 — 100 dyne cm^^ should 
be typical. 



In Table [T] we show the values of T, computed from 
equation ([8|) for different materials assuming Pvap = 1 
dyne cm~^. One can see that these values are consid- 
erably higher than the conventional estimate T, ^ 1500 
K often used in modeling debris disk SEDs. The expla- 
nation for this puzzling difference lies in the fac t that 
the ca nonical estimate is based on the work of .LoddersI 
(|2003f l which explicitly assumes a Solar composition gas 
in equilibrium with sublimating particles to compute Ts . 
Even though a total pressure of 100 dyne cm~^ is as- 
sumed in that work the vapor pressure Pvap of high-Z 
species is going to be much lower in protoplanetary disks 
because of the low abundance of such elements compared 
to H, which contributes m ost to the total pressure. For 
example, in lLodderi (|2003l ) iron was assumed to have an 
abundance (by number, with respect to H) of 3.4 x 10~^, 
which results in vapor pressure of atomic Fe of 3.4 x 10""^ 
dyne cm~^ (assuming that molecular H has dissociated 
and the total pressure in the gas is 100 dyne cm~^). Us- 
ing equation ([8]) we then find Tg « 1300 K instead of 
1600 K typical for a debris disk around a WD, if the lat- 
ter were composed of pure Fe and had a total pressure 
(equal to the Fe vapor pressure) of 1 dyne cm~^. A sim- 
ilar or even larger difference with the canonical estimate 
of Ts arises for other elements listed in Table [TJ 

A deficiency of volatile components (H and He) and 
high relative abundance of metals (can easily be as high 
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as unity) in the gaseous phase of the debris disks around 
WDs naturally results in high values of the (quasi-static) 
sublimation temperature Tg in these systems. This goes 
in the direction of alleviating the puzzle of high temper- 
atures of solid particles inferred from the SED model- 
ing for these objects. However, it does not fully resolve 
this problem because in the simple model of sublima- 
tion Ts is still reached only in the inner rim of the disk, 
i.e. Ts = 7thin(''rim)- The temperature in the bulk of 
the disk just behind the rim is again much lower than 
Ts-. for = 10'' K and Ts = 2100 K as typical for 
olivines (see Table [ij one finds using equation ^ that 
7thick(''rim) ~ 1100 K, which is clearly not enough to 
reproduce the short-wavelength portion of the observed 
SED in many WD-|-disk systems, see Table [2] 

In the following section we provide a complete solution 
to this puzzle by considering particle sublimation in more 
detail. 

3. INNER EDGE STRUCTURE IN THE 
OPTICALLY THICK DISK 

We now present a simple physical model for the struc- 
ture of the inner rim of the disk, the inner part of which 
is optically thick, see Figure [TJ the vertical optical depth 
of such a disk 
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(a and p are the particle radius and bulk density, E is 
the surface density of the disk of solids) is larger than 
unity. This implies that the disk absorbs all incident 
stellar ra.diation because its optical depth to starlight 
( |Rafik2 v|[20 ni) n = t/C ^ 1. According to iRafikovl 
(|2011bD and iMetzger et all (|2012D such optically thick 
inner regions are natural for massive debris disks in which 
aerodynamic coupling to the surrounding gaseous disk is 
strong enough to drive runaway accretion of metals onto 
the WD. In this case the particle surface density at Rm 
is high enough for the inner disk to stay optically thick. 
The key ingredients of our model are 

• the inward migration of particles across the rim re- 
gion, 

• the shielding from starlight provided by the directly 
exposed rim particles to particles further out, 

• the dynamic regime of particle sublimation, as op- 
posed to the quasi-static situation explored previ- 
ously. 

Our primary goal here is to determine the temperature 
inside the rim Trim and the distance Ri„ at which the 
disk gets truncated by sublimation. 

In the previous section, we assumed that particle sub- 
limation occurs at a single temperature Tg , so that parti- 
cles cannot exist in solid form at T > Ts. This, however, 
is not true on time intervals shorter than the time it takes 
to completely sublimate a particle. Using equation (HJ 
one can estimate the instantaneous sublimation timescale 
as the time it takes to completely sublimate a particle 
(neglecting condensation) at a given temperature 
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Figure 1. Schematic representation of the inner rim structure 
and surface density distribution in its vicinity in the optically thick 
(''' t) case. The shaded part of the disk heated to Ti„ < Ts 
receives starlight only at the surface, at grazing incidence angle 
Particles in the inner, unshaded part are directly illuminated 
by the star and are heated to Trim ^5 Ts . Radial width L of this 
exposed rim is determined by equation II16I I. 



Sublimation should be considered as a dynamic (as op- 
posed to quasi-static, like in the previous section) pro- 
cess whenever the particle temperature changes on a time 
scale < ts{Ts). 

When the temperature of a solid object is close to 
its sublimation temperature TsiPvup) (for a given vapor 
pressure Pvap, which at Tg should be equal to P^a*) equa- 
tions ([5|), ([H and (fT3|) allow us to estimate 
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(for a spherical particle of initial radius oq) where oq.i — 
ao/(l cm). Pi = Pvap/(1 dyne cm~^), pi = (1 g cm"^), 
(a)o.i = (q;)/0.1, and p2& = /^/(28mp), as appropriate 
for Si. Because of the rapid scaling of m{T) with T it is 
obvious that ts{T) is a very sensitive function of T, and 
ts{T) < ts{Ts) eve n if T is just slightly higher than Tg. 

As mentioned in ijl.li in the optically thick disk parti- 
cles just outside the rim are shielded from direct starlight 
by the rim particles and their temperature is given by 
Tin — Tthick(-Rm), see equation ^ and Figure [TJ This 
is the highest temperature that one would infer from fit- 
ting the fiat optically-thick disk model to the SED. These 
shielded particles are cool enough for sublimation not to 
be important — an assumption that we check later. 

Particles in the disk migrate inwards due to PR drag 
or aerodynamic coupling to the gaseous disk — this mi- 
gration is what ultimately gives rise to metal accretion 
onto the WD. As particles enter the rim and get exposed 
to direct starlight their temperature rapidly goes up to 
Trim = Tthin(-Rm) > Ti„ (the amouut of energy required 
to heat the particle by several hundred K is small com- 
pared to the heat of sublimation) . According to equation 



5 



([1} the inner rim of the optically thick disk lies at 

i?th'^k^^f^]'. (15) 



The fact that particles in the rim are illuminated by 
virtually unattenuated stellar radiation implies that the 
optical depth of the rim to starlight in the radial direc- 
tion 



0.06 
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(16) 



0.05 



where L is the radial extent of the rim, x is the radial 
distance away from the inner edge of the rim (i.e. the 
location where all particles sublimate; rim corresponds 
to < a; < L), n{x) is the volume number density of 
particles, and we assume the particle cross section for 
starlight to be equal to the geometric cross section tto^ 
(assuming spherical particles). 

Since solids lose mass to sublimation while drifting 
through the rim, particle radius a is a function of x; in 
particular a[x = 0) = 0. The evolution of particle size 
due to sublimation is described by the following simple 
equation 



47r c 



-47ra^m(rrim), 
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which is a simplified version of equation Q in which 
condensation has been neglected. This is a reasonable 
assumption because we will find later that the rim tem- 
perature Trim is significantly higher than the quasi-static 
sublimation temperature Ts{Pva.p)- In this case the flux 
of molecules (or atoms) leaving the particle surface is 
much higher than the flux of particles arriving at it (for 
a given surrounding vapor pressure Pva.p), so that con- 
densation can be neglected. 

We assume that the disk outside the rim is composed 
of particles of a single size oq so that a{x = L) = ao. 
Introducing Vr = dr/dt — dx/dt one can write da/dt — 
Vrda/dx, so that equation (|17|) reduces to 



da _ m{T^ini) 
dx pvr{a,x) 



(18) 



In general Vr{a, x) is a function of both a and x, see e.g. 
equation (I5T]) for the case of PR drag-driven accretion. 

We will now assume that as particles pass through the 
rim and sublimate, their number Qmx Fn does not change 
(until they fully sublimate) even though their mass flux 
varies because their sizes go down as a result of sublima- 
tion. This assumption amounts to neglecting the possi- 
bility of particle breaking or merging during their travel 
through the rim. 

Introducing the vertical thickness of the disk h{x) one 
can use the constancy of Fn to express volume number 
density of particles in the rim n{x) as 



n{x) 



F, 



N 
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We can now express Vr from equation (TTSl) . plug it into 
equation p9]) and substitute the resulting expression for 



C 

Figure 2. Solution of equation l|28|l for ttie scaled temperature of 
the inner rim Trim/T p as a function of the dimensionless parameter 
C given by equation II29I I which contains all information about the 
parameters of the system and particle properties. 

n{x) into the condition (|T6)) to find that 



FnP 
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da 



i) J h{x) dx 
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(20) 



To proceed further we need to make explicit assump- 
tions regarding the behavior of h{x). Debris disks around 
WDs are expected to be similar in properties to dense 
planetary rings around Saturn. The latter have vertical 
thickness comparable to the particle size, which is estab- 
lished by collisions between particles. Thus, it may be 
natural to assume that h(x) ^ a(x), which upon plug- 
ging into equation ((20| and integrating with the condi- 
tion a{L) = ao gives 



FNpal 



4i?„m(T,.im) 



1. 



(21) 



The mass accretion rate of metals onto the WD Mz is 
related to F^ via Mz = F^ x (47r/3)paQ, so that equation 
((2T|) ultimately yields 



M, 



8 27ri?i„ao 



(22) 



One can try another simple approximation for the be- 
havior of h{x), namely assuming that h ao — const. In 
this case one again recovers condition ((22)) with a factor 
of 1/4 instead of 3/8. This similarity of results suggests 
that, for any reasonable assumption regarding the be- 
havior of h{x), the condition 



(23) 



with C ^ 0.05 — 0.1 must be satisfied in the rim. 

Equation (|23|) is the condition that determines the 
value of the inner rim temperature Trim (or, equivalently, 
the inner radius Rin) once the explicit form of m(Ti.im) 
is specified. Given that 2'KRinaQ is the area of the inner 
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rim as seen from the WD, equation (^5)) suggests a sim- 
ple physical interpretation: the disk is truncated at the 
distance Rm , where the full rate of sublimation from the 
area of its inner rim facing the star 27ri?i„ao) roughly 

matches the metal accretion rate through the disk Mz- 
By taking da ao, c?x ~ L in equation we es- 
timate the time tcross it takes particles to cross the rim 
(and sublimate): icross ~ L/vr ~ pao / rniTj-i^) . Using 
equation (|23|) to express rn{Tnui) via Mz and equation 
(fT5|) for Rin we obtain 



tc 



2CMz \Tr 
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where i?+,-2 = i?*/lO^^i?0, Co.i = C/0-1- Interestingly, 
this estimate is independent of the nature of the physical 
process driving particle migration. 

We now plug Rin expressed in terms of Tj-im via equa- 
tion (fT5|) into equation ([23l) to find the following tran- 
scendental equation for Trim only: 



2C 



from which we find 



Tn 



Ari, 



:(lnAi.im) , 

{a)KaRi,aQ 
2CMz 



rim 



(25) 



(26) 



(27) 



with Ai-ini ^ 1. This result can be re- written in the 
following simple form amenable for iterative solution: 



To 



21n 



C = ln 



To 

{a)KQRi,ao 
KMz 



(28) 
(29) 



! 18.6 + In 



(a)o.i-R*,-2ao,iT^,4 
Co.iMz,8 



where Ti,.4 — Ti^jilQ'^ K) and the numerical estimate in 
equation (|29p is done for olivine {Kq = 1.6 x 10^ g~^ 
cm-2 s-\ To = 68, 100 K, see Table p. Figure [D shows 
the exact solution of equation (^5]) for Tri,„/To as a func- 
tion of C . This curve is independent of the system pa- 
rameters and particle properties {Mz, oq, composition), 
which are all absorbed into the definition of C. 

Using equations dS]), (ITT|) and it can be trivially 
shown that 



Ah, 



tc. 



(30) 



As a result, when the time <cross it takes for a particle to 
cross the rim is shorter than the sublimation timescale 
one finds that Ag » Arim and Trim ^ Tg . This illustrates 
our expectation that in the case of dynamical sublimation 
the temperature of particles can be higher than the quasi- 
static sublimation temperature Tg given by equation (|8]). 



?(Ri)- 



R,=Ri„+Ar 



Ar~R,„ 



Figure 3. Schematic representation of the inner rim structure 
and the surface density distribution in its vicinity in the optically 
thin (t|| < 1) inner disk. The inner optically thin part of the disk 
(unshaded) is directly illuminated by the star over a broad range 
of radii Ar ~ Rin- The optically thick part (shaded) starts at 
R\ = Rin + Ar. Compare with Figure [Tl 



For example, for the fiducial values of parameters 
adopted in equation (|29l) one finds for olivine C — 18.6 
and Trim ~ 0.04To « 2700 K, while according to Table 
[T] olivine has Tg « 2100. The inner edge of the disk 
in this case is located very close to the WD surface, at 
Rin w 7i?^, see equation ((T5|) . 

At the same time, just behind the rim the disk temper- 
ature is T,n = Tthick(i?^„) « 1600 K < Tg for = 10^ K, 
see equation ([3]). This verifies our previous assumption 
of relatively low particle temperature (i.e. Tin < Tg) just 
behind the rim, justifying the disregard of particle sub- 
limation in this region. On the other hand, this value of 
Tin is clearly high enough for the disk to produce enough 
near-IR emission corresponding to T ~ 1500 K, in agree- 
ment with observations. 

Equation ((30|) also emphasizes the necessity of particle 
accretion for maintaining the superheated inner rim: if 
Mz — >■ then according to equation ([24)1 icross ^ oo and 
the hot inner rim does not exist. 

4. SUBLIMATION RADII IN DISKS WITH 
OPTICALLY THIN {t\\ < 1) INNER REGIONS 

We now look at the case of a disk, the inner part 
of which is optically thin for incident stellar radia- 
tio n, i.e. ry < 1, see Figure [3l As demonstrated 
bv lBochkarev fc Rafikovl ()2011f ) such situation naturally 
arises for a low mass disk, which starts with T|| < 1 ev- 
erywhere, or for a moderately massive disk, which hasn't 
gone through the runaway accretion phase. In the latter 
case, as shown in iBochkarev fc Rafikovl (|2011t ) , an opti- 
cally thin tail of solid material with Ty ~ 1 (or r ~ ^ ^ 1) 
naturally develops as an inward extension of the optically 
thick part of the disk under the action of the PR drag. 

Particles in such optically thin tail are directly exposed 
to starlight, meaning that their eq uilibriurn temper ature 
is given by equation ([T]). Also, I Rafikovl (|2011bl ) and 
iMetzger et al.l (|2012D have shown that the dynamics of 
these optically thin regions, including the radial drift of 
particles, is determined primarily by PR drag. Then the 
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radial migration speed is just 

3 



1 



Vr.PR 



Stt paoc^ r 



(31) 



The characteristic timescale tpR = r/vr^pR on which the 
particle distance and temperature vary under the action 
of the PR drag is then 



Stt pa^c o 



10 yr- 



0.2i?f: 



(32) 



Since tppj is much longer than the sublimation 
timescale ts given by equation (fT3|) it is clear that in 
the optically thin disks sublimation must be occurring in 
a quasi-static fashion: particles slowly drift inward under 
the action of the PR drag and their temperature steadily 
rises. At some radius -R*Jj'" their temperature reaches Tg, 
and particles turn into metal gas on a (short) sublimation 
timescale tg . That means that the inner edge of the opti- 
cally thin disk is set by the condition Tthin(^*Ji'") = Ts, 
with Ts given by equation ([SJ. Thus, 



nthin _ ( ^ 



(33) 



In particular, according to Table [T] we need to take Ts ~ 
2100 K for olivine, which when plugged in the equation 
((33)) yields w lli?^ for = lO'' K. This is about 
60% further from the star that in the case of an optically 
thick disk, see ^ 

5. APPLICATION TO OBSERVED SYSTEMS. 

We now apply ideas developed in fJSl to a sample of 
observed WDs with debris disks. We start by rewriting 
the expression (^5]) for C as 

(34) 
(35) 



C — -|- Cp 



C.=ln 



R*T^ 



M7 



Cp = In 



{a)KQaQ 



2CTi 



Here C^, is a parameter, which depends only on mea- 
surable properties of the system — WD radius, effective 
temperature, and metal accretion rate. All parameters 
characterizing the particle properties — i^o, ao, etc. — 
are absorbed into Cp. Assuming a particular composi- 
tion of particles and a value of particle radius oq fixes Cp 
and allows one to obtain a theoretical relation between 
Trim and C* using equations (|26l) . (p4| . and ([35]). By 
looking at different particle compositions one can com- 
pare the corresponding theoretical T'rim(C*) curves with 
the properties of observed systems. 



Such comparison requires the knowledge of i?*, T^,, 
Mz, which we take from the literature. One also needs to 
know Trim for each of the WD-|-disk systems, and we de- 
rive this parameter as Ttiiin from equation ([S]), in which 
we use Tin — the innermost disk temperature inferred 
from the SED fitting — for Tthick- We use the val- 
ues of Tin determined in the literature when available, 
and we provide our own fits otherwise. The summary of 
WD-|-disk parameters used in our comparison with the- 
ory is provided in Table [2l 



In Figure|3]we show theoretical Trim(C'*) curves for dif- 
ferent particle compositions. In our calculations we al- 
ways assume ao = 1 cm particles, (a) = 0.1, and C, = 0.1 
(all dimensional quantities are expressed in CGS units). 
We also plot the locations of observed systems from Ta- 
ble [2] in Ci, — Trim space with hexagons. 

As expected, very refractory particles made of 
graphite, SiC, and AI2O3 are characterized by consid- 
erably higher values of theoretical Trim (for the same C*) 
than if they were to have more volatile compositions, e.g. 
were made of iron. The difference in Trim can easily ex- 
ceed 10^ K. 

The vast majority of observed systems lies in between 
the two extremes determined by the iron and graphite. It 
is clear from this plot that the pure graphite composition 
is not acceptable for particles in the observed WD+disk 
systems — all of them are below the corresponding the- 
oretical curve. Also, only a handful of systems lie close 
to the theoretical Trim(C'*) curve for iron. The majority 
of observed WD-|-disk systems tend to gravitate towards 
Trim(C'*) curves computed for CAI and olivine-like com- 
positions. At the same time about a third of the systems 
in the upper right corner of the figure are consistent with 
more refractory compositions such as SiC or AI2O3. 

When comparing characteristics of observed systems 
with theoretical predictions for Trim(C'*), a couple of is- 
sues have to be kept in mind. First, observational deter- 
mination of parameters of the WD-|-disk systems is prone 
to errors. This is not so serious for the determination of 
Ti,, which is typically quite accurate, or i?^,, which does 
not span a large range anyway. However, the determi- 
nation of Mz from the data depends on the unknown 
composition of the parent body that formed the disk, 
and may have large error bars. On the other hand, 
depends on these characteristics only logarithmically, so 
that even large uncertainties in these parameters would 
result in a relatively small horizontal shift of observa- 
tional points in Figured 

The uncertainty in measuring Trim is much more seri- 
ous. This is because the determination of Tin relies on 
fitting the fiat disk model to the SED, and Tin can be 
highly dege nerate with other p arameters, such as the disk 
inchnation (|Girven et al.ll2012l ). Also, according to equa- 
tion dSJ Trim oc Tf,/^, so that the errors in determination 
of Tin from SED directly propagate into the uncertainty 
in Trim. As a result, observational data points in Figure 
|3]can have significant vertical errorbars. 

Another thing to keep in mind, is that when computing 
the theoretical Trim(C'^) curves we make certain assump- 
tions about particle properties, such as their size ao or 
accommodation coefficient {a). Variation of these pa- 
rameters from their adopted values affects the value of 
Cp and causes horizontal shift of the Trim(C^,) curves. 
For example, increasing the value of accommodation co- 
efficient (a) from 0.1 to 1 displaces the theoretical curves 
to the left by AC* = 2.3. This would put observational 
data points in better agreement with the more refractory 
particle compositions. 

6. DISCUSSION. 

The physical model for the inner rim structure in the 
optically thick case presented in fj3] naturally allows us 
to explain the high inner disk temperatures T^ inferred 
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Table 2 

Properties of disk-hosting WDs used in this work 



Name 


SpT 


M^, 






logio M 


Gas Disk 


T 


T ■ 

rim 


C* 


Ref. 






Mq 


Rq 


K 


g/s 


Detected 


K 


K 


CGS 




GD 16 


DAZB 


0.59 


0.014 


11500 


8.0 




1300'' 


2460 


21.0 


1 


GD 133 


DAZ 


0.59 


0.014 


12200 


8.5 




1200 


2380 


19.9 


1.2 


GD 40 


DBZ 


0.59 


0.013 


15200 


9.9 




1200 


2560 


17.1 


1 


GD 56 


DAZ 


0.60 


0.015 


14200 


8.5 




1700 


3160 


20.3 


1,2 


J1228-I-1040 


DAZ 


0.77 


0.011 


22020 


9.3 


yes 


1670 


3610 


19.0 


3,4 


PGlUlo-f lol 


DAZ 


0.61 


0.014 


19300 


9.3 




1200 


2770 


19.0 


1,2 


Ton345 


DBZ 


0.70 


0.010 


18600 


9.4 


yes 


1500 


3180 


18.4 


5,6,7 


SDSS1043+0855 


DAZ 


0.66 


0.012 


17900 


9.0 


yes 


1400 


3000 


19.4 


8 


G29-38 


DAZ 


0.62 


0.013 


11700 


8.7 




1200 


2350 


19.3 


1,9,10 


GD 362 


DAZB 


0.73 


0.013 


10500 


10.4 




1200 


2260 


15.2 


1,9,11 


SDSS0959 


DAZ 


0.64 


0.012 


13280 


7.9 




1600 


2970 


21.3 


12 


SDSS1221 


DAZ 


0.73 


0.011 


12250 


7.7 




1400 


2640 


21.6 


12 


SDSS1557 


DAZ 


0.42 


0.018 


22810 


8.8 




1400 


3250 


20.8 


12 


GD 61 


DBZ 


0.71 


0.011 


17280 


8.81 




1300 


2820 


19.7 


13,14 


J0738-I-1835 


DBZ 


0.84 


0.010 


13950 


11.11 


yes 


1600 


3020 


13.9 


16 


HE 0110-5630 


DBAZ 


0.71 


0.012 


19200 


8.4 




1000 


2450 


20.9 


17,18 


HE 1349-2305 


DBAZ 


0.67 


0.012 


18200 


8.7 




1700 


3430 


20.1 


17,18 



^ These numbers were calculated in this paper. 

b (1) [F arihi ct al. ^OOl; (2) Jura ct al. 2007a; (3) Brinkworth ct al. 200 91: (4') IGansicke et al][200a: (5) 
Farihi et al. 2 010: (6) McUs ct al. 2010; (7) Gansickc ct al. 2008; (8) Brinkw orth et al.| 20ia (9)IFari hi et al ' 



20081: (10)T2uckcrman ct al. 2003; (11) Zuckcrman ct al. 2007; (12) Farihi ct al. 2012; (13) Farihi et ai 



201TI: (14) [Jura &: Xu.2012.: (15) iKilic et al..,2012, : (16) .Dufour et al..i2012i ; (17) iGirven et al.., 20121 : (18) 
Koester et al.ll200a . 



t=0.1 
<a> = 0.1 
an= 1 cm 



3000 



"^2500 

.^iC o 

: Fe 

2000 - AlgOg 

graphite 

- olivine 




1500 



10 



15 

dn(R.T?/MJ 



Figure 4. Comparison of observed WD+disk properties with 
theo retical predictions in C*— T^im space (C* is defined by equation 
H35| lV Theoretical Tj-\^{Ci,) curves computed for different particle 
compositions are labeled on the plot. Their calculation assumes 
ao = 1 cm particles, f = 0.1, and (a) = 0.1. Note that most of 
the observed systems are consistent with particles being made of 
Si-bearing materials, such as olivine or CAI. 



from the SEDs of debris disks around some WDs. The 
existence of a narrow inner rim of the disk heated to 
a temperature Trim above the quasi-static subhmation 
temperature Tg (see equation (|H])) is the key ingredient 
of the model. 

The radial width of the inner rim L can be estimated 
by multiplying the time to cross it tcross by the velocity 
Vr.pR, given by equations (1241) and (j31|) correspondingly: 



10 cm 



ao,i-R*,-2i*, 



Z,8 



0.2AU^ 



Ri. 



(36) 



Thus, one typically finds the width of the inner rim to 
be ~ 10 particle radii. 

Note that the radial speed of particles in massive disks 
can be affected by aerodynamic coupling between the 
particulate and gaseous disks, and in consequence devi- 
ate from Vr.pR. Nevertheless, equation (j36l) serves as a 
reasonable order of magnitude estimate of L and clearly 
demonstrates that L ^ As a result, the contribu- 

tion of the hot inner rim to the SED of the debris disk 
is completely negligible, and its spectrum is determined 
only by emission from the parts of the disk located be- 
hind the rim. 

Our results in |j3][4]allow us to address the differences in 
spectra of disks with optically thick or thin inner regions. 
In Figure [5] we show several spectra produced by disks 
around a T, = 10^ K, = O.Oli?© WD located 10 pc 
away from us, and inclined with respect to our line of 
sight with cosi = 0.5. The model, which is optically 
thick everywhere, has constant optical depth r = 10 and 
extends from the outer radius Rout = ^0 to R*^'^^ « IRi, 
given by equation ([T5|) . Models with optically thin tails 
also have constant optical depth r = 10 between Rout = 
Rq and some intermediate radius which is different 
for each model. Inside of Ri we assume an optically thin 
tail with constant r — ({Ri) (or T|| = r/Ri) to extend 
from Ri down to i?*^'" « Hi?* given by equation ((33)) . 
This is the characteristic distribution of r in the inner 
optically thin tails of th e disks e volving under th e action 
of the PR drag, see iBochkarev fc Rafikovl ()2011[ ). 

One can see that the disk which is optically thick ev- 
erywhere produces more flux. This is expected because 
disks with optically thin tails do not extend as far in- 
ward, and are inefficient at absorbing and re-radiating in 
regions interior to r = Ri. The spectral shape is also 
different, in part because particles in the optically thin 
tail are hotter than particles in the optically thick tail 
at the same radius. This may allow one to diagnose the 
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Figure 5. Spectra of debris disks with optically thick and opti- 
cally thin inner regions. All models feature optically thick regions 
with T = 10, which extend from Rout = Rq to in the opti- 

cally thick case and to r = Ri (indicated in the panel) in the mod- 
els with optically thin tails. Inside Ri tails have r = C(^l) 1- 
Stellar parameters are indicated on the panel. See text for more 
details. 

presence of an inner optically thin tail using just the disk 
SED. Such optically thin tails may be expected in sys- 
tems characterized by Mz ~ 10^ g s~^, i.e. close to the 
value pr ovided by PR drag a lone. In the runaway sce- 
nario of iMetzger et al.l ()2012[ ). one expects systems with 

higher Mz to be evolving due to aerodynamic coupling 
with the gaseous disk, in which case the disk is optically 
thick all the way down to i?-^"^"^- 

Comparison of our theory with characteristics of ob- 
served WD-|-debris disk systems shows that in general 
(barring the uncertainties related to measurement errors 
and poorly constrained modeling parameters) properties 
of these systems are consistent with Si-rich particle com- 
position. In other words, we find that CAI- or olivine- 
like compositions of particles are in reasonable agreement 
with the locations of the inner rims in the majority of ob- 
served disk-hosting systems. 

This result reinforces previous conclusions about the 
Si-rich nature of the accreted material based on different 
and independent lines of evidence. In particular (and 
most importantly), direct measurements of the metal 
abundances in the WD atmospheres show that the com- 
position of accreted material is consistent with that of 
the inner Solar Sy stem bodies , which are known to be 
Si-rich ( Zuckerman et al.l 120071 : iKfein et al.l[20Tol 2011; 
iJuraeLaDiQl^ These measurements also demonstrate 
the accreted bodies to be carbon-poor purall2006[ ). which 
is again consistent with our results — essentially none of 
the observed WD+disk systems lie close to the C-based 
curve in Figure HI Additional evidence in favor of Si-rich 
particle composition comes from the measurement of 10- 
fxm bump in debris di sk spectra obtained with Spitzer 
IRS pura et al.ll2"007H 2009). This feature is usually in- 
terpreted as being produced by the /im-size silicate par- 
ticles. 




0.03 cm 



30 cm 



10 15 20 25 

C. = ln(R.T?/M2) 

Figure 6. Effect of varying the particle radius a on the theoretical 
Trim—C* dependence for olivine and comparison with observational 
data. 

Using these independent lines of evidence supporting 
the Si-rich nature of the debris disk constituents we may 
approach our findings from a different perspective. In 
particular, by postulating disk particles to be Si-rich we 
can put constraint on their sizes. Results presented in 
Figure |4] do a reasonably good job at reproducing char- 
acteristics of observed systems by assuming a — I cm 
particles. Varying a would displace the theoretical curves 
horizontally and they would remain consistent with ob- 
servations only within a certain range of particle sizes. 
Figure [6] illustrates this variation of the Trim — C** re- 
lation; one can easily infer from it that particle sizes 
should lie within the range a = 0.03 — 30 cm. Other- 
wise the properties of the inner disk rims in the majority 
of the observed systems will not be consistent with our 
theoretical calculations. 

Interestingly, this range of particle sizes is consistent 
with other indirect measu rements of a reported in the 
literature. In particular, iGraham et al.l (jlQQOf ) found 
a < 10 cm based on the varia bility of the repro cessed 
IR emission of disk particles. Me tzger et al.l (|2012f ) found 
a < several cm to provide the best fit to the runaway pic- 
ture of the disk evolution. Finally, Saturn rings, which 
are thought to be rather close in properties to circum- 
WD debris disks, ar e also predominan tly composed of 
1 - 100 cm particles (|Cuzzi et al.ll2010D . 

Our model naturally explains the presence of solid par- 
ticles even around hot WDs, with Ti, w 20, 000 K, e.g. 
J1228+1040 and SDSS1557. Conventional theory finds 
it difficult to account for such systems. Indeed, equation 
^ predicts that around = 20, 000 K, = O.OlSi?© 
WD directly illuminated particles must have a temper- 
ature of 1700 K at the tidal radius of ^ Rq. This is 
significantly higher than the sublimation temperature of 
1300 — 1500 K usually assu med based on protoplanetary 
disk studies (|Loddersll2003[ ). Our calculations first show 
that in fact the sublimation temperature Tg can easily 
be higher than 1700 K, see Table [TJ which guarantees 
the survival of even the optically thin disks with directly 
exposed particles within tidal radii of hot WDs. Second, 
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in the optically thick case, shielding of the disk by the in- 
ner rim particles allows Ri„ to be as small as O.3i?0 (for 
Ri, = 10~^i?o, and keeping all other parameters equal 
to their values in equation (HH)). 

Our present calculations were designed to demonstrate 
the main qualitative features of the inner rim structure 
and thus made a number of simplifying assumptions. 
One of them is the single size of particles in the disk, 
while in reality a distribution of particle sizes should be 
present. We expect that in this case the value of a in the 
definition pSI) of Cp would be replaced with some prop- 
erly weighted average of the particle size distribution, 
but the main results would not change. 

Another simplification is the assumed single chemical 
composition of all particles. If the disk contains particles 
of different compositions, with different Kq and Tq, then 
one may expect a "multi-rim" structure to form, in which 
different chemical species sublimate at different radii. In 
this case the inner radius of the disk would be determined 
by the properties of the most refractory particles in the 
disk that survive at the closest separation from the WD. 
Observations of the inner disk properties (i.e. Tin) would 
then be sensitive to characteristics of only this particular 
particle population (as long as the disk is optically thick 
everywhere) . 

7. SUMMARY. 

We explored the structure of the inner parts of com- 
pact debris disks around WDs with the goal of resolv- 
ing the "inner rim puzzle" — the difficulty with rec- 
onciling the high inner disk temperatures inferred from 
the SED with the material properties of putative con- 
stituent particles. We first show that because of the 
much higher vapor pressure of metals in these hydrogen- 
poor disks compared to the hydrogen-rich protoplanetary 
disks, the quasi-static sublimation temperature Tg of dif- 
ferent species in circum-WD disks is typically 300-400 
K higher than in their conventional protoplanetary ana- 
logues. This revised value of determines the (smaller 
than was thought before) value of the inner radius for 
the optically thin disks, given by equations ([8]) & (|33|) . 

We demonstrate that optically thick circum-WD disks 
feature narrow inner rims, which are superheated above 
Ts- This allows inner disk radii in such systems, de- 
scribed by equations (fTSi) . p6)) . & p7|. to lie quite close 
to the WD, easily at separations ^ lOi?^,. The main 
physical ingredients needed for the existence of such su- 
perheated inner rim are (1) accretion of particles through 
the disk, which can be easily maintained at the necessary 
level by Poynting-Robertson drag, (2) shielding of parti- 
cles behind the rim from starlight by the rim particles, 
and (3) dynamic nature of the sublimation process in- 
side the hot rim. The combination of these ingredients 
naturally allows particles to reach temperatures of or- 
der 1600 — 1700 K just behind the rim, which is needed 
to explain the SEDs of some systems. Particles inside 
the rim are heated to 2500 — 3500 K and undergo rapid 
sublimation as they migrate in. Using this model we 
can naturally explain the existence of particulate debris 
disks even around hot WDs, with effective temperature 
> 20,000 K. 

We compare our predictions with existing observations 
of the WD-|-disk systems. We find that properties of par- 
ticles in debris disks are consistent with Si-rich compo- 



sition, such as olivine or CAI-like material. Very refrac- 
tory (such as graphite) or more volatile (such as iron) 
compositions are clearly disfavored by this comparison. 
Assuming that circum-WD disks are indeed composed 
of Si-rich particles we constrain typical particle size to 
lie roughly between 0.03 and 30 cm, in agreement with 
other indirect evidence for cm-size objects in such disks. 
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APPENDIX 

THERMAL BALANCE IN THE DISK OF SOLIDS. 

Grains in the inner rim are being heated and cooled by four processes: (1) heating by starhght, (2) heating by gas, (2) 
coohng by thermal radiation from particle surfaces, and (4) removal of thermal energy by sublimating atoms/molecules. 
All these processes scale linearly with the surface area of the particles. As a result, the temperature of grains directly 
exposed to starlight (assuming full absor ption of the inc oming radiation) is implicitly given as a function of the distance 
from the WD by the following formula ()Podolakll20Tol ): 

+Qgas = 'TT'i + gsub, (Al) 

Qsub = m{T)Lsuh, Qgas = ^^2^^^^ '^'^ga^ ~ ' 

where Lgnh is the specific heat of sublimation of the particle material, s is the efficiency of heat exchange between gas 
and particles, and pg is the gas density. 
Using equation ([TU)) we can estimate Sg « Mz / {Sniy) so that 

M . iMl, (A3, 

This allows us to compare the contribution of gas heating Qgas with stellar irradiation (first term in the left hand 
side of equation (jAip ) : 

Qgas 2s GM^Mz ( T \ eMz,sM.,i ( ^ 

Therefore, gas heating is typically unimporta nt for t he thermal balance of particles compared to heating by starlight, 
in contrast to the conclusion reached by Jur a et al.l ([2007b), who looked at conduction in gas phase as the means to 
lower particle temperatur e. This difference i s predominantly caused by the high gas density (^ 10^ times higher than 

in equation (jASp ) used in lJura et al.l (|2007bl ). 

Using prescription (O with (a) = 0.1 and the typical (for olivine) value of Lgub — 3.2 x 10^° erg g^^ from lKimura et all 
Hoo^ we can also estimate the relative contribution of sublimation to the cooling of particles by evaluating Qsub/cT . 
We find this ratio to be about unity for particles heated to « 3400 K. At T = 3000 K the ratio of the energy loss by 
sublimation to radiative cooling rate is about 0.1. Thus, for WD+disk systems with Trim ^ 3000 K one can safely 
neglect Qsub in equation (jAll) . Then the thermal balance everywhere in the disk is determined by the equilibrium 
between stellar heating and radiative cooling only, which provides justification for using equations ([T]) and ^ in this 
work. This assumption is good for the majority of observed systems shown in Figure 21 and even for a handful of 
systems with Trim ~ 3000 — 3500 our theoretical curves should still be at least qualitatively correct. 

DATA ON THE MASS SUBLIMATION RATES 

Here we present the details on the derivation of mass sublimation rates for different elements shown in Table [1] 
Olivines Calculation of the vapor pressure for the olivine-like silicate grains (e.g. Mg2Si04) is compli cated due to 
the fact that these molecules do not exist in the gas phase. Nevertheless. iGuhathakurta &i Draind (|1989( ) suggest the 
following expression for the (number) rate of Si sublimation from the olivine surface: i?si « 7 x 10^°(a) exp(— 68, 100/T) 
cm~^ s~^. The mass sublimation rate of olivine is then given by HonRsu where = 141 rup is the mean molecular 
weight of Mg2Si04. 

Graphite For pure graphite IGuhathakurta fc Draind (|1989( ) give the (number) rate of C sublimation from the 
graphite surface of Rc ~ 4.6x 10'^°(a) exp(— 81, 200/T), which then allows us to calculate Kq from the mass sublimation 
rate -R qMCi where iig = 12 mp is the mean molecular weight of carbon. 

CAI iRichter et al.l (|2007D consider evaporation of CAI-like liquids and come up with the following (number) rate 
of Si escaping a CAI-like surface: i?si ~ 2.5 x lO^^(a) exp(-69, 400/T). Using gehlenite (Ca2Al2Si07) as a typical 
CAI-like material (mean molecular weight 274 nip) we obtain sublimation parameters indicated in Table [TJ 

Iron iZatisev et al.l (|200in provide the data on the vapor pressure of Fe: P^^^ = 2.8 x 10" exp(-45, 400/T) Pa 
for T 1800 — 1900 K. From these data we determine the mass sublimation rate according to the formula rfiFe — 
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(a)P^p (/iFc/27rfcBr)^^^, where /^fg ~ 56mp and we take T = 1600 K with the expectation that the thermophysical 
parameters of Fe remain roughly the same at this t emperature. 

SiC Using the thermophysical data presented in iChase et al.l (|1985f) we derive the following fit to the behavior of 
the vapor pressure of Si above the SiC surface in the range T = 1800 - 3000 K: Pfl^ = 9 x lO" exp(-73, TOO/T) 
dyne cm~^. Since the surface loses C atoms in addition to Si we evaluate SiC mass loss rate as wsic = 
(a)P4'p (Aisic/Msi) (Msi/2^fcBT)'/', where /igic = 40mp, ^gi - 28mp and we take T = 2400 K. 

AI2O3 For pure corundum (AI2O3) the basic reaction which is thermodynani ically most likely is AI2O3 — » 2A10+0 
(a different reaction dominates for AI-AI2O3 mixture, see lBrewer fc Searcvlll951[ ). Using the data in lChase et al.l (|19850 
we find that the vapor pressure of O above the corundum surface is = 3.6 x 10^^ exp(— 80, 500/T) dyne cm~^ for 
T = 1800 — 3000 K. Accounting for the mass of Al leaving the surface together with O we arrive at the sublimation 
characteristics of corundum presented in Table [TJ 



